Astronomy & Astrophysics manuscript no. paper 
December 4, 2012 



©ESO 2012 



Letter to the Editor 



What controls the magnetic geometry of M dwarfs? 

T. Gastine^, J. Morin^, L. Duarte^, A. Reiners^, U. R. Christensen^, and J. Wicht^ 



' Max Planck Institut fiir Sonnensystemforschung, Max-Planck-StraBe 2, 37191 Katlenburg Lindau, Germany 

e-mail: gastine@mps.mpg.de 
^ Institut fiir Astrophysik, Georg-August-Universitat Gottingen, Friedrich-Hund Platz, 37077 Gottingen, Germany 



>: 

^ ■ 

o: 

^ ■ 
(N. 



X: 



December 4, 2012, Revision : 1.107 



(N 
O 

o 



o 



ABSTRACT 



Context. Observations of rapidly rotating M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole- 
dominated and multipolar geometries. In dynamo models, the relative importance of inertia in the force balance - quantified by 
the local Rossby number - is known to have a strong impact on the magnetic field geometry. 

Aims. We aim to assess the relevance of the local Rossby number in controlling the large-scale magnetic field geometry of M dwarfs. 
Methods. We explore the similarities between anelastic dynamo models in spherical shells and observations of active M-dwarfs, 
focusing on field geometries derived from spectropolarimetric studies. To do so, we construct observation-based quantities aimed to 
reflect the diagnostic parameters employed in numerical models. 

Results. The transition between dipole-dominated and multipolar large-scale fields in early to mid M dwarfs is tentatively attributed 
to a Rossby number threshold. We interpret late M dwarfs magnetism to result from a dynamo bistability occurring at low Rossby 
number. By analogy with numerical models, we expect different amplitudes of difi^erential rotation on the two dynamo branches. 
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1. Introduction 

M dwarfs - the lowest-mass stars of the main sequence - are 
of prime interest to study stellar dynamos operating in phys- 
ical conditions quite remote from the solar case. During the 
past few years, their surface magnetic fields have been inves- 
tigated using two complementary approaches (for recent re- 
views see Donati & Landstreet 2009; Reiners 2012). On the one 
hand, with spectroscopy in unpolarized light (i.e. measuring the 
Stokes parameter 1) the average magnetic field strengths (here- 
after referred to as (Bj)) of tens of stars spanning the whole 
M-dwarf spectral class have been derived. From such measure- 
ments Reiners & Basri (2007) have concluded that no signifi- 
cant change in (Bj) occurs when stars become fully convective 
(spectral type M3/M4) and Reiners et al. (2009) could constrain 
the rotation-magnetic field relation ((Bj) increases towards short 
rotation periods until it reaches a plateau at ~ 3 kG). On the 
other hand, the measurement of circular polarization (Stokes V) 
in spectral lines of M dwarfs provides a constraint on the orien- 
tation and polarity of the field at large and intermediate scales. 
Modelling of the surface field performed with Zeeman-Doppler 
imaging (ZDl, Semel 1989) on a sample of about 20 M0-M8 
dwarfs points towards a broad variety of magnetic field geome- 
tries. Partly-convective stars as well as a few fully-convective 
ones feature complex magnetic structures (e.g. Donati et al. 
2008; Morin et al. 2010) while most fully-convective ones host 
a strong axial dipole component (Morin et al. 2008a,b). 

Explaining such a diversity in the magnetic field geometry 
is one of the main goals of stellar dynamo theory. However, 
as the numerical simulations cannot directly access the parame- 
ter regime where such dynamos are thought to operate, asymp- 
totic scaling laws are of prime interest to check the relevance 
of the dynamo models (Christensen 2010). Scaling laws de- 
rived from geodynamo-like models (i.e. incompressible flow and 



no-slip boundaries) successfully fit the magnetic field strength 
of a broad range of astrophysical objects, encompassing Earth, 
Jupiter and some rapidly-rotating M dwarfs (Christensen et al. 
2009). Such similarities between the dynamos in planets and 
rapidly -rotating stars motivate the comparison of some planetary 
dynamo results with observations of active M dwarfs. 

In geodynamo models, the relative contribution of inertia and 
Coriolis force in the global balance is known to have a strong im- 
pact on the magnetic field geometry (e.g. Sreenivasan & Jones 
2006). Christensen & Aubert (2006) suggest that the ratio of 
these two forces can be quantified by the so-called "local Rossby 
number" defined by Ro/ - Wims/fi/, / being the typical flow 
lengthscale. This dimensionless parameter has been found to 
be a rather universal quantity that allows to separate dipolar 
and multipolar dynamo models: a sharp transition between these 
two types of dynamo occurs around Ro/ ^0.1. However, this 
has been recently challenged by some models that employed 
stress-free mechanical boundary conditions, more appropriate 
when modelling stellar dynamos (e.g. Goudard & Dormy 2008). 
Simitev & Busse (2009) for instance found that a dipolar and 
a multipolar magnetic field can coexist at the same parameter 
regime depending on the initial condition. This magnetic bista- 
bility phenomenon can lead to multipolar solutions even for 
Ro, < 0.1 (Schrinner et al. 2012). 

In addition, most of these studies have been conducted un- 
der the Boussinesq approximation that assumes the reference 
state to be constant with radius. While this assumption is suit- 
able for modelling the geodynamo, it becomes questionable in 
the stellar interiors, where density increases by several orders 
of magnitude. Anelastic and compressible dynamo models in- 
deed indicate a strong influence of the density contrast on the 
geometry of the magnetic field: some weakly stratified simu- 
lations of Dobler et al. (2006) are dipole-dominated, while the 
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Strongly stratified models of Browning (2008) are multipolar. 
The parameter study of Gastine et al. (2012) extends the bistabil- 
ity phenomenon to a broader parameter range and confirms that 
strongly stratified models tend to produce multipolar dynamos. 

Following the conclusions by Schrinner et al. (2012) and 
Gastine et al. (2012) (hereafter GDW12), the main aim of this 
letter is to compare the results of the later study with spectropo- 
larimetric observations of M dwarfs. 

2. The dynamo model 

We consider MHD simulations of a conducting anelastic fluid in 
spherical shells rotating at a constant rotation rate Q. about the 
z-axis. Following Gilman & Glatzmaier (1981), the governing 
MHD equations are non-dimensionalised using the shell thick- 
ness d - To - ri as the reference lengthscale and Qr^ as the time 
unit. Our dynamo model results are then characterised by sev- 
eral dimensionless diagnostic parameters. The rms flow veloc- 
ity for instance is given by the Rossby number Ro = Ui-^s/^d. 
Following Christensen & Aubert (2006), we also employ the 
aforementioned local Rossby number Ro/ = Roiu/n, that is 
known to be a more appropriate measure to quantify the impact 
of inertia on the magnetic field geometry. The mean spherical 
harmonic degree 4 is obtained from the kinetic energy spectrum 
and relates to the typical flow lengthscale / through: 

l^nd/l, with 4^y/ ^"''"f , (1) 

<«•«) 

where ue is the flow at a given spherical harmonic degree ( and 
the brackets correspond to an average over time and radius. 

The magnetic field strength is measured by the Elsasser num- 
ber A - B^^^l PIJ.AD., where p is the density, and and A are 
the magnetic permeability and difFusivity. The geometry of the 
surface field is quantified by its dipolarity /dip = B^^j /,,=o('" ~ 

'"o)/ Tj^pm ^?m('" ~ '"o)' '^he ratio of the magnetic energy of the 
dipole to the magnetic energy contained in spherical harmonic 
degrees up to ^njax = 11 (see Christensen & Aubert 2006). 

The dimensionless MHD equations are advanced in time 
with the spectral code MagIC (Wicht 2002; Gastine & Wicht 
2012) that uses the anelastic formulation of Lantz & Fan (1999) 
and has been validated against several dynamo benchmarks 
(Jones et al. 201 1). We rely in the following on the results of the 
parameter study of GDW12. Although these simulations have 
been initially tailored to model the dynamo of giant planets, the 
differences in the reference properties (i.e. density and gravity 
profiles) are not expected to be of prime influence. Similarities 
between dynamos in planets and low-mass stars are indeed em- 
phasised by geodynamo-based scaling laws (Christensen 2010) 
and previous Boussinesq models employed in the stellar dynamo 
context (e.g. Schrinner et al. 2012; Simitev & Busse 2012). 

3. Spectropolarimetric observations 

Spectropolarimetric observations of 23 active M0-M8 dwarfs 
with rotation periods ranging from 0.4 to 19 days have been car- 
ried out. For each star at least one time-series of unpolarized and 
circularly polarized spectra sampling a few rotation periods has 
been obtained. The data reduction and analysis is detailed by 
Donati et al. (2006, 2008) and Morin et al. (2008a,b, 2010). 

The relative importance of inertia with respect to the Coriolis 
force in the convection zone of these stars is assessed through an 



empirical Rossby number given by 

ROemp = , (2) 

^conv 

where Tconv is the empirical turnover timescale of convection 
based on the rotation-activity relation (Kiraga & Stepien 2007). 
This Rossby number misses explicitly the flow lengthscale / in- 
volved in Ro/. However, as Tconv is based on the average convec- 
tive turnover time it encompasses this scale information to some 
extent. We thus use Roemp as our best available proxy for Ro/. 

For each obtained spectrum, an average line profile with in- 
creased signal-to-noise ratio is computed using the least-squares 
deconvolution technique (LSD, Donati et al. 1997). Each time- 
series of LSD profiles is modelled with ZDI, resulting in a map 
of the large-scale component of the surface magnetic field vector 
that satisfies a maximum-entropy criterion. The large-scale mag- 
netic fields of most of these stars fall into two distinct groups: 
one is dominated by a strong axial dipolar component and the 
other by a much weaker and non-axisymmetric field. We could 
not identify any ZDI reconstruction bias that would spuriously 
produce such a separation, and exclude an effect of limited reso- 
lution as both groups span largely overlapping ranges of v sin /. 

Following Morin et al. (201 1), we define an Elsasser number 
based on the averaged unsigned large-scale magnetic field (By) 
which roughly characterises the ratio of Lorentz and Coriolis 
forces, rj is obtained by rescaling a reference solar magnetic dif- 
fusivity of 770 = 10" cm^s"' (Stix 1989) to the M dwarfs with 
a MLT argument. We also consider the fraction of the magnetic 
energy that is recovered in the axial dipole mode in ZDI maps. 
The spatial resolution of such maps mostly depends on v sin / and 
the actual degree and order ^^ax up to which the reconstruction 
can be performed ranges from 4 to 10, although very little en- 
ergy is recovered in modes with A < C < 10. We therefore di- 
rectly compare this quantity to the dipolarity employed in nu- 
merical models and term them both /dip. We however note that 
in simulations, /dip does not strongly depend on the chosen ^niax, 
whereas for the observation-based dipolarity, considering the ra- 
tio of magnetic energy in the axial dipole relative to the total 
magnetic energy derived from unpolarized spectroscopy (instead 
of the large-scale magnetic energy derived from spectropolari- 
metric data with ZDI) would lead to much lower values of /dip 
(cf. Reiners & Basri 2009). We attribute this diff'erence to the 
low magnetic Reynolds number (Rm ~ 100 - 500) accessible 
by numerical simulations which does not allow for a significant 
small-scale field to be generated - hence the weak dependence of 
/dip on /"max - while in stellar interiors large-scale and small-scale 
dynamo action likely coexist (e.g. Cattaneo & Hughes 2009). 

4. Results 

Figure 1 shows /dip versus Ro; in the numerical models, which 
use the anelastic approximation with density contrasts up to 
Np - ln(pbot/ptop) - 3 and stress-free (or mixed) boundary con- 
ditions. Figure 2 displays the relative dipole strength of M stars 
against Rogmp derived from spectropolarimetric observations. 

The numerical models cluster in two distinct dynamo 
branches: the upper branch corresponds to the dipole-dominated 
regime (/dip > 0.6), while the dynamo models belonging to 
the lower branch have a multipolar magnetic field (/dip < 0.2). 
In agreement with the geodynamo models, the dipolar branch 
is bounded by a maximum Ro/, beyond which all the mod- 
els become multipolar (Christensen & Aubert 2006). Although a 
small dependence on the shell thickness is noticeable, the dipo- 
lar branch vanishes around Ro/ ^0.1 for both spherical shell 
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Fig. 1. /dip plotted against Ro/ in anelastic dynamo models. Red (grey) 
symbols correspond to simulations in thick (thin) shells (r,7ro = 0.2 and 
r-Jra = 0.6, respectively) and their size is scaled according to the value 
of the surface field, expressed in units of the square root of the Elsasser 
number. Each type of symbols corresponds to a density contrast. The 
two closed symbols correspond to two cases further discussed in Fig. 3. 
The vertical lines mark the tentative limits for dipolar dynamos. 
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Fig. 2. Observational counterpart of Fig. 1 . The vertical dashed line 
marks the tentative limit for the dipolar regime. For the two stars ex- 
hibiting the largest temporal variations, the individual epochs are shown 
and connected by a vertical red line. Dotted red circles with black ar- 
rows correspond to stars from Morin et al. (2010) for which a definite 
ZDI reconstruction could not be achieved, but an upper limit for the ro- 
tation period and an estimate of dipolarity were derived. The two closed 
symbols correspond to two selected stars discussed in Fig. 4. 



geometries considered by GDW12. However, in contrast with 
most of the previous Boussinesq studies, the multipolar branch 
also extends well below Ro/ ^0.1, where both dipolar and mul- 
tipolar solutions can coexist (see also Schrinner et al. 2012). In 
the parameter range covered by our numerical models, bistability 
of the magnetic field is in fact quite common, meaning that both 
a dipole-dominated and a multipolar field are possible configura- 
tions at the same set of parameters. The type of solution is then 
selected by the geometry and the amplitude of the initial seed 
magnetic field (Busse & Simitev 2006; Simitev & Busse 2009). 
The multipolar branch found at Ro/ < 0.1 is partly composed 
by the multipolar component of such bistable dynamos but also 
encompasses all the stratified models with Np > 2, for which no 
dipolar solutions are found (GDW12). 

The comparison between these results and the observations 
of stellar magnetic fields suffers from the difficulty to relate the 
diagnostic parameters of the dynamo models (i.e. Ro/, A and 
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Fig. 3. Surface radial magnetic field 6,-(r = tq) and axisymmetric zonal 
flows for a dipolar dynamo model with Np = 1.7 (a) and a multipolar 
one with Np = 2 (b). The maps of 6,. have been low-pass filtered up to 
^max = 10. Magnetic fields are expressed in units of the square root of 
the Elsasser number and velocities in units of the Rossby number. 



/dip) to their stellar counterparts. Within these limits, the separa- 
tion into two dynamo branches seems to be relevant to the sam- 
ple of active M dwarfs displayed in Fig. 2. In fact, all the early M 
stars (with > 0.5 Mq) show multipolar magnetic fields with 
/dip < 0.2; being slow rotators they fall into the Roemp > 0.1 
regime (Reiners & Mohanty 2012). The observations of mid M 
dwarfs suggest a possible transition between dipole-dominated 
and multipolar magnetic fields close to Rogmp ~ 0.1, although 
CE Boo does not fit into this picture (green square in the upper 
right corner of Fig. 2). Late M dwarfs (with < 0.15 Mq) 
seem to operate in two different dynamo regimes: the first ones 
show a strong dipolar magnetic field, while others present a 
weaker multipolar magnetic structure with a significant time 
variability (emphasised by the vertical red lines in Fig. 2). These 
important time variations of the dipole strength are also fre- 
quently observed in multipolar dynamo models with low Ro/ 
(e.g. Schrinner et al. 2012, GDW12) but are not visible in Fig. 1, 
where time-averaged properties are considered. We note that 
the values of Tconv entering Roemp are poorly constrained for M 
dwarfs. Assuming a stronger (weaker) mass-dependence would 
mostly expand (shrink) the x-axis of Fig. 2 without changing our 
conclusions. The uncertainties on /dip typically lie in the range 
0.1-0.3 (see discussion in Morin et al. 2010), which does not 
question the identification of two distinct branches. 

The two dynamo branches found in our numerical models 
also differ by their main force balance, at least in the bistability 
region (i.e. Ro/ < 0.1). The dipolar branch encompasses mod- 
els with Elsasser number around unity that suggest a first-order 
contribution of the Lorentz force in the main balance (the dy- 
namo then operates under the so-called "magnetostrophic bal- 
ance", e.g. Sreenivasan & Jones 2006). In contrast, at low Ro/, 
the models belonging to the multipolar branch have weaker mag- 
netic fields (see Fig. 1 and Tab. 2 in GDW12), meaning that the 
Lorentz force may not enter in the first-order balance. In that 
case, significant geostrophic zonal flows (i.e. aligned on coax- 
ial cylinders) can possibly develop (e.g. Gastine & Wicht 2012). 
Note that beyond Ro; > 0.1, the multipolar dynamos are strong 
enough to yield a magnetostrophic force balance due to larger 
Rayleigh numbers. This correspondence between high dipolar 
fraction and strong large-scale magnetic field is also present in 
spectropolarimetric observations (Morin et al. 201 1). Figure 3 il- 
lustrates these two types of dynamos for two typical models with 
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Fig. 4. Surface radial magnetic field of V374 Peg (a) and GJ 1245 B 
(b) recovered with ZDI from spectropolarimetric observations. The 
field has been reconstructed up to fy„.,^ = 10 (4) for V374 Peg 
(GJ 1245 B). Surface dilferential rotation of V374 Peg has been derived 
by Morin et al. (2008b) from spectropolarimetric observations, while 
this was not possible for GJ 1245 B. Magnetic fields ai'e expressed in 
units of the square root of the Elsasser number. 

Ro/ < 0. 1 . While the upper panel shows a dipole-dominated so- 
lution with a very weak differential rotation, the lower one has a 
multipolar field that goes along with significant zonal flows that 
become nearly geostrophic close to the equator These strong 
zonal winds affect the dynamo mechanism via the f2-effect, i.e. 
the production of toroidal field by the shear. This Q-effect plays 
only a minor role in the field production of the dipole-dominated 
models that can be categorised as "o'^ -dynamos", according to 
the mean-field description (e.g. Chabrier & Kiiker 2006). In con- 
trast, dynamos on the multipolar branch can be classified as aO. 
or a^Q, at least in the low Ro/ regime. 

To a certain extent, these differences are also noticeable in 
the reconstructed fields of M dwarfs as illustrated in Fig. 4. 
While V374 Peg (upper panel) has a dipole-dominated magnetic 
field and a very weak differential rotation, GJ 1245 B (lower 
panel) presents a weaker amplitude multipolar field with impor- 
tant time variability (one of the vertical red lines in Fig. 2). If, as 
suggested in the present study, the multipolar fields observed in 
several late M dwarfs are the consequence of a dynamo bistabil- 
ity occurring at low Ro/, the numerical models would then sug- 
gest a significant differential rotation in these stars (AQ/fi ~ 5% 
in the Fig. 3b model). However, this cannot be verified from 
the available data as surface differential rotation has been deter- 
mined for one star of the sample only (a dipolar one, see Fig. 4a). 

5. Conclusion 

Spectropolarimetric observations of active M dwarfs and dy- 
namo models show a broad variety of magnetic geometries (see 
Gastine et al. 2012; Morin et al. 2010, and references therein). 
In both cases, dipolar and multipolar large-scale magnetic fields 
are found to coexist at low Rossby numbers. In the present letter 
we critically discuss the analogy between these two results. 

We derive observation-based quantities aimed to reflect the 
diagnostic parameters employed in the numerical models (Ro;, 
A and /dip), although these crude proxies are not expected to 
provide a direct quantitative match. Within these limits, we draw 
an interesting analogy between the observational parameters and 
their numerical counterparts: for large values of the Rossby num- 
bers multipolar fields are found, while below a critical value 



around Ro/ ~ Roemp ~ 0.1, a bistable region exists where both 
dipolar and multipolar fields can be generated. 

Several limitations must be noted though. ( i) The spectropo- 
larimetric sample is biased as all stars at high (low) Roemp are 
partly (fully) convective. Thus it is not yet clear if the change in 
/dip observed around Roemp ~ 0.1 can be attributed to a thresh- 
old in Ro/ or rather to the drastic changes in stellar structure 
occurring at the fully-convective transition. ( ii) As the numerical 
models of GDW12 do not attempt to model a tachocline, they 
might miss some important features of early M dwarfs mag- 
netism. However these issues do not question the validity of 
the agreement between observations and simulations regarding 
the existence of a bistable dynamo regime at low Ro/ for fully- 
convective stars. (Hi) In numerical models, the dipolar branch 
only exists for mild density contrasts (Np < 2), much below 
the stratification of stellar interiors. Different assumptions from 
those used by GDW12 could possibly extend the dipolar regime 
towards higher stratifications, for instance by using different 
values of Prandtl numbers (Simitev & Busse 2009) or radius- 
dependent properties (e.g. thermal and ohmic diffusivities). 

However, the analogy between dynamo simulations and 
magnetic properties of M dwarfs can be further assessed with 
more realistic numerical models, and additional observations as 
it implies that: ( i) stars with multipolar fields can be found over 
the whole parameter range where also dipole-dominated large- 
scale fields are observed; (ii) in the bistable domain, stars on the 
multipolar branch have a much stronger surface differential ro- 
tation than stars hosting dipole-dominated large-scale fields. 

Acknowledgements. TG and LD are supported by the Special Priority Program 
1488 (PlanetMag) of the German Science Foundation. JM is funded by a post- 
doctoral fellowship of the Alexander von Humboldt foundation. AR acknowl- 
edges research funding from DFG grant RE 1664/9-1. 

References 

Browning, M. K. 2008, ApJ, 676, 1262 

Busse, F H. & Simitev, R. D. 2006, GAFD, 100, 341 

Cattaneo, F & Hughes, D. W. 2009, MNRAS, 395, L48 

Chabrier, G. & Kiiker, M. 2006, A&A, 446, 1027 

Chiistensen, U. R. 2010, Space Sci. Rev, 152, 565 

Christensen, U. R. & Aubert, J. 2006, GJl, 166, 97 

Christensen, U. R., Holzwarth, V., & Reiners, A. 2009, Nature, 457, 167 

Dobler, W., Stix, M., & Brandenburg, A. 2006, ApJ, 638, 336 

Donati, J.-F, Forveille, T., Cameron, A. C, et al. 2006, Science, 311, 633 

Donati, J.-F & Landstreet, J. D. 2009, ARA&A, 47, 333 

Donati, J.-F, Morin, J., Petit, P, et al. 2008, MNRAS, 390, 545 

Donati, J.-F, Semel, M., Carter, B. D., et al. 1997, MNRAS, 291, 658 

Gastine, T., Duarte, L., & Wicht, J. 2012, A&A, 546, A19, [GDW12] 

Gastine, T. & Wicht, J. 2012, Icarus, 219, 428 

Gilman, R A. & Glatzmaier, G. A. 1981, ApJS, 45, 335 

Goudard, L. & Dormy, E. 2008, Europhysics Letters, 83, 5900! 

Jones, C. A., Boronski, P, Brun, A. S., et al. 201 1, Icarus, 216, 120 

Kiraga, M. & Stepien, K. 2007, Acta Astronomica, 57, 149 

Lantz, S. R. & Fan, Y. 1999, ApJS, 121, 247 

Morin, J., Donati, J., Petit, P, et al. 2008a, MNRAS, 390, 567 

Morin, J., Donati, J.-F, Forveille, T., et al. 2008b, MNRAS, 384, 77 

Morin, J., Donati, J.-F, Petit, P, et al. 2010, MNRAS, 407, 2269 

Morin, J., Dormy, E., Schiinner, M., & Donati, J.-F 2011, MNRAS, 418, L133 

Reiners, A. 2012, Living Reviews in Solar Physics, 9, 1 

Reiners, A. & Basri, G. 2007, ApJ, 656, 1121 

Reiners, A. & Basri, G. 2009, A&A, 496, 787 

Reiners, A., Basri, G., & Browning, M. 2009, ApJ, 692, 538 

Reiners, A. & Mohanty, S. 2012, ApJ, 746, 43 

Schrinner, M., Petitdemange, L., & Dormy, E. 2012, ApJ, 752, 121 

Semel, M. 1989, A&A, 225, 456 

Simitev, R. D. & Busse, F H. 2009, Europhysics Letters, 85, 19001 

Simitev, R. D. & Busse, F H. 2012, ApJ, 749, 9 

Sreenivasan, B. & Jones, C. A. 2006, GJI, 164, 467 

Stix, M. 1989, The Sun. an Introduction (Berlin: Springer- Verlag) 

Wicht, J. 2002, Physics of the Earth and Planetary Interiors, 132, 281 



4 



